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An assessment of the experimental uncertainty in obtaining the kinetic activation energy from ther- 
mogravimetric analysis (TGA) data is presented. A neat phenolic resin, Borden SC 1008, was heated at 
three heating rates to obtain weight loss vs temperature data. Activation energy was calculated by two 
methods: the traditional Flynn and Wall method based on the slope of log(^) versus 1 JT, and a modifi- 
cation of this method where the ordinate and abscissa are reversed in the linear regression. The modified 
method produced a more accurate curve fit of the data, was more sensitive to data nonlinearity, and gave 
a value of activation energy 75% greater than the original method. An uncertainty analysis using the 
modified method yielded a 60% uncertainty in the average activation energy. Based on this result, the 
activation energy for a carbon -phenolic material was doubled and used to calculate the ablation rate in 
a typical solid rocket environment. Doubling the activation energy increased surface recession by 3%. 
Current TGA data reduction techniques that use the traditional Flynn and Wall approach to calculate 
activation energy should be changed to the modified method. 


Nomenclature 

A = pre-exponential factor in Arrhenius equation, Eq. (6) 
B r = bias uncertainty of an experiment result 
C = degree of conversion or weight loss, Eq. (3) 

E - kinetic activation energy, Eq. (6) 

F - integral weight-loss function, Eq. (10) 

/ = conversion weight function, Eq. (8) 

k - decomposition reaction rate constant, Eq. (6) 

N = number of results or data points 
n = decomposition reaction order, Eq. (8) 

P r = precision uncertainty of the experimental result 

q = heating rate 

R = universal gas constant 

r = linear regression correlation coefficient 

S = precision index 

T = absolute temperature 

t = time 

U r - uncertainty in the experimental result 
W = sample weight at time t 
X = mean result of a sample set 
X, = single result within a sample set 
T = volume fraction of resin 

p = material density 
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Subscripts 

CAL = calibration value 

CF = curve-fit value 

j, k - resin components 

/ = fibrous component 

m = material component 

o = original or initial value 

r = residual or final value 


Introduction 

P OLYMERIC and polymer-based composite materials are 
widely used in the aerospace industry for bonding, insu- 
lating, and sealing. A common application of these materials 
is in ablative cooling of combustion chambers and rocket noz- 
zles. In this application a matrix of high-temperature-resistant 
fibers (such as carbon or silica) is impregnated with an organic 
phenolic or epoxy resin. The impregnated matrix is cure hard- 
ened into an appropriately shaped rigid structure. When ex- 
posed to the high-temperature characteristics of the burning 
propellant gases, the material undergoes a complex process 
combining chemical decomposition, melting, sublimation, 
charring, and evaporation. 1 3 This physicochemical process of 
ablation causes the formation of a char layer on the surface, 
which is a poor heat conductor, thus forming a protective ther- 
mal blanket for the underlying structure. 2,4 

The Space Shuttle solid rocket motor (SRM) nozzle thermal 
protection system is composed of an ablative carbon— phenolic 
composite material. Current material research seeks to reduce 
the amount of composite required in the nozzle and thereby 
reduce nozzle weight (and thus launch costs). Conservative, 
i.e., thicker, ablative liners are often used because of the degree 
of uncertainty currently present in ablation simulation codes. 
Improved analytical techniques for predicting material ablation 
are thus required. 56 

Thermogravi metric analysis (TGA) 7 is a technique particu- 
larly suited for quickly characterizing the thermal degradation 
of a material with time and/or temperature. Kinetic reaction 
parameters as well as weight fractions of material components 
are calculated from the TGA data and input to computational 
ablation simulation codes. As with any type of experimental 
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approach (and resulting calculations) some error is inherent in 
the process. Assessing the experimental uncertainty associated 
with the TGA process would lead to increased confidence in 
the resultant simulation-code predictions of nozzle ablation. 

Many studies have been conducted that compare various 
methods of determining kinetic reaction parameters of poly- 
meric materials from TGA data. 8-12 These methods include 
curve fitting (using empirical TGA data and iteratively solv- 
ing for reaction parameters), approximating the integral of 
the TGA decomposition equation, and obtaining kinetic pa- 
rameters from integral weight fractions as a function of tem- 
perature. Although no consensus has been reached as to the 
most accurate method, 13 the procedure developed by Flynn 
and Wall 8,14 is widely used and has been incorporated into the 
data-reduction programs of several commercially available 
TGA instruments. This method uses linear regression analysis 
to determine the relationship between log heating rate and re- 
ciprocal temperature. The slope of the plot resulting from the 
regression equation is used to calculate activation energy. 

Typically, kinetic reaction data are generated at low heating 
rates. Because extremely high heating rates are observed in 
actual rocket nozzles the generation of TGA data at these high 
rates are also being examined. 15 The specialized equipment that 
has been developed for these studies is not commercially avail- 
able. A comprehensive literature search has indicated that no 
One has systematically assessed the experimental uncertainty 
associated with the determination of kinetic reaction param- 
eters by one particular methodology and commercially avail- 
able instrument. In the study presented in this paper typical 
procedures for determining activation energy from TGA data 
were evaluated. The uncertainty analysis methods of Coleman 
and Steele 16 were then employed to perform a detailed analysis 
of the propagation of uncertainty through the experimental de- 
termination of activation energy. The value of these methods for 
determining uncertainties associated with linear regression anal- 
ysis has been demonstrated using Monte Carlo simulations. 17 

Most data for carbon-phenolic composites are restricted by 
the International Traffic in Arms Regulations. 18 Because the 
purpose of this research was to assess the uncertainty in the 
experimental TGA methodology by examining typical proce- 
dures for determining activation energy, neat (no fibers) poly- 
meric resin SC 1008 was used. Without the effect of carbon 
fibers it is not expected that the activation energy of a neat 
phenolic resin would be the same as that of a phenolic com- 
posite. The uncertainty in activation energy of a composite 
may be greater than that of a neat resin because of errors 
caused by sample nonuniformity. Ten trials of the neat poly- 
meric resin were conducted at each of three different heating 
rates. Based on the resultant calculated experimental uncer- 
tainty, standard activation energies for a carbon-phenolic ma- 
terial were increased by an equivalent factor and used to cal- 
culate the ablation rate in a typical solid rocket environment. 

Ablation Process 

Several computer codes have been written to predict the 
ablation rate of carbon-phenolic insulation materials in a 
rocket nozzle environment. The charring material thermal re- 
sponse and ablation program (CMA) code, and derivatives 
thereof, is currently used by NASA and its contractors for 
thermal analysis of the SRM ablative nozzle. The physico- 
chemical ablation model is shown schematically in Fig. 1. 
Heating causes chemical decomposition in upper layers of the 
composite material, which results in the formation of a char 
layer at the top of a pyrolysis zone. Volatile chemical species 
released in the pyrolysis zone percolate through the char to the 
surface. The only effect on the virgin material is a small loss 
of volatiles. The backup or substrate material is unchanged by 
the effects of heat from the ablation process. 

The ablation mechanism is assumed to be a surface chemical 
equilibrium, diffusion-controlled phenomena, where external 
combustion products react with the carbonaceous char and py- 


> FREE 

Chemical STREAM 

BOUNDARY 
LAYER 

CHAR 


PYROLYSIS 
ZONE 

VIRGIN 
MATERIAL 

BACKUP 
MATERIAL 

Fig. 1 Schematic of material ablation process model. 

roly sis gases at an exponential rate. 12 The ablation rate is pro- 
portional to the diffusion rate and concentration of reacting 
species that diffuse through the boundary layer. For pyrolysis 
of a composite material the resin filler is assumed to consist 
of two or more components that decompose separately (des- 
ignated by j and k ), whereas the reinforcing fibrous material 
(designated by l) is the third component that can decompose. 
The instantaneous density of the composite is thus given by 4 

P = T( ft + p k ) + (1 - I>, (1) 

The volume fraction of the resin T must be input into the 
simulation codes. It is assumed that each of the three com- 
ponents decompose according to the kinetic reaction described 
by 



The decomposition rate constant k m must be accurately deter- 
mined to calculate the material ablation rate. 

Thermogravimetric Analysis 

The material ablation process is quantified in terms of the 
conversion weight fraction C, which is the fraction of material 
already converted or decomposed at any instant in time. 9,19 For 
a given sample of material, the conversion weight fraction is 
defined by 



C=(W 0 - W)/{W a - W r ) (3) 

The decomposition rate of a material is a function of a tem- 
perature-dependent rate constant and the degree of conversion 
(or weight loss) 8,19 


— = k(T)f(C) (4) 

df 

The temperature is assumed to be a linear function of time and 
the specified heating rate 

T = qt + const (5) 

The rate constant is assumed to be represented by the Arrhe- 
nius form 8,14,20 


k(T) = A exp(—E/RT) (6) 


Substituting Eqs. (5) and (6) into Eq. (4) yields the funda- 
mental kinetic model for the thermogravimetric decomposition 

, 8 20 21 
rate 


dC 
d T 


= -/(C)exp 
<7 



( 7 ) 
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The TGA system monitors the weight of a sample as it is 
heated isothermally or at a specified heating rate. The isother- 
mal approach is often avoided because it is impossible to in- 
stantaneously raise the material to the desired temperature be- 
fore some initial weight loss occurs. 7 Therefore, the specified 
heating rate approach is generally used. The most common 
output from the TGA system is a thermogram, which is a plot 
of the percent of material remaining versus temperature. Figure 
2 shows a typical thermogram for phenolic resin SC 1008 
heated at a rate of 20°C/min. Each region of slope change 
(indicated as A, B, and C) on the curve corresponds to an area 
of rapid weight loss. Typically, these areas represent different 


reactions U king place in the material as it is heated. By taking 
the time dt rivative of the curve, as plotted in Fig. 3, the dif- 
ferent reac ions can be seen more clearly as inverted peaks. 

If the reiction rate continues to decrease with decomposi- 
tion, the nte of conversion is proportional to the amount of 
remaining londegraded material, 9 14 19 described by 

/(C) = (1 - C) n (8) 

Because reaction parameters are usually obtained with com- 
mercial instruments from the beginning of the decomposition. 
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a first-order reaction model is commonly assumed (n = l), 22 
and the decomposition rate may be expressed as 


d C 
d T 



- C)exp 



(9) 


A quick, direct method for determining activation energy 
from TGA data was proposed by Flynn and Wall. 14 They as- 
sumed that A , /(C), and E are independent of 7, and that A 
and E are also independent of C. As discussed in Refs. 14 and 
21, separating variables and integrating Eq. (9) yields 



where x = EfRT. The integral in Eq. (10) cannot be evaluated 
directly and the value of p(x) must be estimated. Taking the 
log of Eq. (10) yields 

log[F(C)] = lo g(AE/R) - log (q) + log p(x) (11) 

Doyle 21 found that for E/RT ^ 20, log p(x) may be approxi- 
mated at a constant percent conversion by 


log Ip(EfRT)] = -2.315 - QA51EIRT (12) 


Doyle’s calculations show that the value of the coefficient, 0.457, 
does not vary more than ±1% over the range 29 ^ E/RT < 
42. Substituting Eq. (12) into Eq. (11) and differentiating at a 
constant percent of conversion, i.e., C = const, yields 



Fig. 4 Schematic of a typical TGA instrument. 1, balance; 2, 
furnace; 3, sample; 4, counterweight; 5, gas inlet; and 6, gas exit. 

Parameters Affecting Experimental Uncertainty 
of the TGA Process 

There are several sources of uncertainty in any TGA anal- 
ysis. Among these are fluctuations in the balance, measurement 
of temperature by the thermocouple, and sample preparation. 
The weight measurements are affected by vibrations, static 
electricity, the thermocouple touching the balance pan, and 
purge-gas flow rate. From experiment trials and other 
studies, 824-26 it was determined that these factors are relatively 
minor and controllable. 

Because TGA samples are measured in milligrams, small 
sample sizes can potentially present a problem in terms of 
sample content homogeneity if proper precautions are not 
taken. With composite material the resin and fiber content may 
vary too widely in small samples. Samples for this investiga- 
tion are pellets taken from a sheet of cured neat resin. Because 
no fibers are present sample homogeneity is less of a concern. 


d(log q) = 


-Q.457F 

R 



Solving Eq. (13) for E yields 


-R d log(g) 
0.457 [ d(l IT) 


(13) 


(14) 


The bracketed term in Eq. (14) can be obtained from TGA 
data as the slope of a plot of log(<?) versus MT at a constant 
degree of conversion. 

Thermogravimetric Measurements 
TGA Instrument 

A Perkin-Elmer TGA-7 instrument 23 was used in this re- 
search to generate thermograms for the neat resin decompo- 
sition. This instrument records the thermal response of a ma- 
terial as it is heated at a specified heating rate. A typical TGA 
is shown schematically in Fig. 4. The TGA-7 primarily con- 
sists of an ultramicrobalance and a furnace element. Using a 
null-balance design the microbalance detects changes in 
weight as small as 0.1 pg. Its maximum capacity is 130 mg. 
According to the manufacturer balance accuracy is better than 
0.1% of full scale. 23 The furnace has the capability to operate 
from ambient temperature to 1000°C at the minimum linear 
heating rate of 0.1°C/min and a maximum heating rate of 
200°C/min (in 0.1°C/min increments). This feature accom- 
modates the requirement of the Flynn and Wall method for a 
minimum of three distinct heating rates to calculate activation 
energy. 14 Other major components of the system include a ther- 
mocouple, gas purge system, and a data-acquisition system. 


Thermocouple Calibration 

A chromel-alumel thermocouple is employed for the TGA 
temperature measurement. Gradients between the sample and 
the thermocouple should be minimized or at least quantified. 
A Curie-point calibration 21,23,2427 was performed for each heat- 
ing rate that was used. This calibration consists of heating from 
2 to 5 ferromagnetic metal standards over a specified heating 
range at the desired heating rate. The standards are lowered 
into the furnace with an appropriately shaped magnet around 
it. Once heated these ferromagnetic metal standards lose their 
magnetism at a specific temperature and this loss in magnetism 
is manifested as an apparent weight loss that is measured by 
the TGA. In 1978 a study was conducted by the International 
Confederation for Thermal Analysis 28 to evaluate the variance 
in standard temperature with instrument and technique. Eight- 
een different researchers using instruments from nine different 
manufacturers participated in the study. The bias in standard 
temperature was found to be markedly dependent upon magnet 
placement, temperature, and thermocouple location in a par- 
ticular instrument. Based on these results TGA instrument 
manufacturers have since sought to minimize the bias as a 
result of this experimental technique with specific instructions 
on the placement of the magnet and thermocouple for their 
systems. 

The temperature distribution within the TGA furnace also 
varies with heating rate. It has been observed that the variance 
associated with heating rate increases near the temperature 
limit of the instrument. Consequently, a calibration should be 
performed for each combination of thermocouple, sample 
holder, and heating rate in accordance with the test objectives. 
It is often recommended to perform the Curie calibration using 
two standards. 23 29 However, four Curie standards were used 
for calibrating the TGA-7 in the research of this paper. Be- 
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cause the TGA-7 instrument software only accepts two cali- 
bration values the two input values were taken from a linear 
regression that was based on the four calibration points. This 
was done to minimize error from a two-point fit and to allow 
calculation of the associated bias from the linear relationship 
between standard and thermocouple temperatures. 

Experimental Procedure 

For this study, 10 pellet samples of neat polymeric resin, 
Borden SC 1008, were heated to 1000°C at each of three dif- 
ferent heating rates of 5, 10, and 20°C/min. These heating rates 
were chosen based on manufacturer’s recommendations. 22 In 
addition, temperature errors associated with the use of pellets 
is minimized with low heating rates. Argon was used as the 
sample purge gas at a flow rate of 20 cc/min. The samples 
were cured at 121°C and ambient pressure. To minimize errors 
caused by sample nonuniformity, variations in resin pellet 
weight were kept as small as possible. The TGA manufacturer 
recommended that the variance in sample weight be no more 
than ±20%. 23 For this study sample weights were held within 
±0.6 of 10 mg or ±6%. To eliminate moisture samples were 
held at 100°C for 10 min and then reweighed. This new dry 
weight was input to the TGA as a basis for the weight-loss 
thermograms. The amount of moisture per sample was less 
than 2% by mass. Data from seven previous trials of the same 
material at 20°C/min 30 were combined with this data set and 
used to calculate activation energy. Temperature data at each 
percent conversion for each thermogram were obtained by us- 
ing TGA-7 data analysis software. The percent conversions 
chosen were based on a manufacturer’s recommendation 22 that 
data be taken from the beginning of the weight- loss curve. The 
material first starts to degrade at this point and practical ma- 
terial lifetime ends. 23 It is this activation energy that is reported 
for this study. In addition, for the SC 1008 resin used, activa- 
tion energy could not be computed beyond 6% conversion. For 
this material, several reactions occur and other techniques 
would have to be employed to determine activation for each 
subsequent reaction. 8,9 


The brack ited term in the denominator is the slope of a plot 
of 1 IT ver ;us log(^). The activation energy was calculated by 
both the standard Flynn and Wall method [Eq. (14)] and mod- 
ified method of Eq. (15). The modified method was used in 
the uncertainty analysis. 

Experimental Uncertainty Analysis 

In this nvestigation the experimental uncertainty analysis 
procedures detailed by Coleman and Steele 1633 are employed. 
In this ap)roach the overall uncertainty in an experimental 
result r is defined in terms of bias and precision uncertainties 
by 

U r = (B] + P 2 ) 1 ' 2 (16) 

The bias i ncertainty in the final result is an estimate of the 
fixed error within the experimental system that cannot be min- 
imized by averaging the results of several trials. Precision un- 
certainty is a measure of the randomness or scatter in the data. 
This type of uncertainty arises as a result of measurement sys- 
tem precis: on errors and also variations in test conditions from 
trial to tria of the same experiment. The true result lies within 
the interva : r ± U r 95% of the time. 

The unc;rtainty in activation energy was calculated based 
on the modified form of Eq. (15). Because Doyle’s coefficient 
of 0.457 is not truly constant, the variable P is substituted in 
its place. Thus, E is a function of four experimental variables 

E = £(/?, P, r, q) (17) 

Following standard TGA procedures 23 the final value of E is 
assumed t< be the average value of activation energy at 1,3, 
and 6% ccnversions, E^ 6 . Based on Eq. (15) the average ac- 
tivation energy is thus a function of the average of the slope 
of the 1/T versus \og(q) curve at 1,3, and 6% conversion 


P( SLOPED + SLOPE,* + SLOPE,*) 


Calculation of Activation Energy from TGA Data 

The variance in sample temperature with heating rate forms 
the basis of the calculation of activation energy by the Flynn 
and Wall technique. This approach, which is commonly used 
in the data reduction software of commercial TGA instruments, 
assumes a linear relationship between the ordinate log q and 
the abscissa 1/T. 14 A linear least-squares method is used to 
determine the slope of this plot. This slope is input to Eq. (14) 
to determine activation energy. 

However, in a typical TGA experiment, temperature is the 
experimental value that is actually being measured as opposed 
to heating rate that is essentially fixed. Therefore, temperature 
is the dependent variable, while heating rate is independent. 
By definition, the objective of the least-squares method is to 
determine the equation of the line that minimizes the vertical 
distance of each data point from that line. 31,32 The ordinate (or 
uncontrolled variable) y in an experiment depends upon the 
abscissa (or controlled variable) x that has minimal variation. 
In the TGA experimental method temperature is obviously the 
uncontrolled variable and heating rate is controlled. Therefore, 
it would seem logical to minimize the variance in temperature 
rather than heating rate. Again, this is the opposite of the ap- 
proach used in the software of many TGA manufacturers. 

Rearranging Eq. (14) such that temperature is the dependent 
variable results in 


There am a total of 1 1 1 different values of 7\ each measured 
at a speed ed corresponding heating rate. Therefore the full 
functional orm for the average activation energy is 

P ■> Pi* * - • T m , <7i, . . . ^m) (19) 

Each expei imental parameter in E has an associated bias rep- 
resented by B r , B P , B t , and B q , respectively. Engineering con- 
stants are lypically known with greater accuracy than experi- 
mental parameters and, thus, the bias in the gas constant B R is 
assumed negligible. Bp is taken from Doyle 21 as ± 1% of 0.457. 
According to Perkin -Elmer B q is negligible. The bias in tem- 
perature, h >wever, is not negligible. The elemental sources of 
bias considered in the temperature measurements are B CF as a 
result of carve fitting the thermocouple calibration data, and 
B std from he Curie standards. Following Ref. 16, B CF is ob- 
tained fron 


&cf — 2 ( 5 cf ) 


where 


S C f — 


N 


2 - w] 


( 20 ) 


( 21 ) 



id(log q\ 


(15) 


The value of B C f was calculated to be 7°C. The bias in the 
standard was estimated to be of the same magnitude and is 
taken to be 7°C. These elemental error sources are combined 
by the rss method to obtain the total temperature bias as 


0.457 


Br = (B 2 C F + B 2 STD ) ] ' 2 = 9.9°C 


( 22 ) 
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It is also necessary to consider the effects of correlated terms 
on uncertainty. These terms arise when elemental errors in dif- 
ferent measured variables originate from identical sources. For 
the TGA measurements of this study the temperature readings 
are correlated because each measurement of temperature is af- 
fected by the same bias uncertainties in the standards used in 
the calibration procedure. These covariance terms are approx- 
imated by the sum of products method.* 4 

Considering the various sources of bias uncertainty relevant 
to this study the resultant bias limit in activation energy is 


Bl = 



B: 


III / „\ 112 / ■ 

( dE\ v' ( dE 

§ 2 (aT* J J?. [dT, 


BtJ 


(23) 


A forward finite differencing method 35 is used to calculate the 
partial derivative terms that are called sensitivity coefficients. 

The precision limit P r specifies the confidence interval 
around a single result, where the biased mean of the experi- 
mental result lies 95% of the time. 33 The precision uncertainty 
includes both correlated and uncorrelated terms. Although 
there has been one documented study that indicated correlated 
precision terms can be significant, 36 these terms are typically 
assumed to be insignificant and are neglected in this paper. 
The precision uncertainty is estimated for large sample sizes 
of 10 or more. Also, the result of activation energy is averaged 
from several trials. Therefore, the precision uncertainty in the 
average activation energy is 33 

Pi = 2 (S e /Vn) (24) 


ward finite differencing technique. From the partial derivatives 
the systematic and correlated bias uncertainties were then cal- 
culated by Eqs. (23-25). The total experimental uncertainty in 
E was then determined by Eq. (26). 

Activation Energy Calculation 

The Flynn and Wall method has been widely used in TGA 
instrumentation to calculate activation energy. This method 
utilizes the relationship between the log of the heating rate and 
reciprocal temperature to derive an equation for activation en- 
ergy. According to Flynn and Wall 14 this relationship can best 
be described by a linear mathematical model. A summary plot 
of log(<?) versus 1 IT data for various percent conversions is 
shown in Fig. 5. A linear regression method was used to de- 
termine the slope at each percent conversion. The sample r is 
a dimensionless quantity that describes the degree of linear 
association between two variables. 35 For the curve-fits in Fig. 
5, r ranges from 0.91 to 0.60 for conversion percentages of 
1 -6%, respectively. The calculated slopes are input to Eq. (14) 
to determine E. 

In a TGA experimental setting the heating rate is fixed and 
the temperature is the observed response that is being mea- 
sured. By definition, the ordinate (uncontrolled) variable in an 
experiment depends upon the abscissa (controlled) variable. 
The objective of the least-squares method is to determine the 
equation of the line that minimizes the vertical distance of each 
data point from the regression line. 31 ' 32 There is a negligible 
variation in heating rate throughout a specified heating profile. 
However, temperature, for a particular percent conversion at a 
fixed heating rate, varies slightly as shown in Fig. 6 for 1% 
conversion. Comparable variations in temperature were ob- 
served for each of the other percent conversions. Because this 
scatter is present in the temperature data, the effect of inter- 


where 

.- N -i 1/2 

*■ = [aTTT 2 ® - *’>J <»> 

Based on the calculated bias and precision uncertainties, the 
uncertainty in the mean activation energy E is 

Ue = (Bl + Pi) 1 ' 2 (26) 

Results and Discussion 

The effects of uncertainty in activation energy determined 
by TGA for the decomposition of a polymeric neat resin was 
assessed. Ten samples of neat resin (Borden SC 1008) were 
heated to 1000°C for each of three different heating rates to 
obtain a thermogram of weight loss versus temperature. Acti- 
vation energies obtained from the Flynn and Wall method and 
a modification of this method, where the ordinate and abscissa 
in the regression are reversed, were each examined. Addition- 
ally, uncertainties associated with the experimental approach 
were investigated. Finally, the effect of uncertainty in activa- 
tion energy on the ablation rate for a composite solid rocket 
nozzle was also investigated. 

A Fortran program was developed to calculate activation 
energy by both the traditional Flynn and Wall and modified 
methods. 37 The program was structured to read temperature 
data for all three heating rates, convert to the appropriate units, 
and store the converted data in an array. The converted data 
were initially used to calculate E for each percent conversion 
based on Eq. (14). Activation energies calculated by the For- 
tran program were compared to the E values calculated by the 
Perkin -Elmer instrument software. These values generally 
agreed to within ± 3 kJ/g-mol or 5%. Therefore, the activation 
energies presented are those calculated by the Fortran program 
and not the Perkin -Elmer software. The Fortran program cal- 
culates the partial derivatives with respect to each temperature 
variable and Doyle’s coefficient defined in Eq. (23) by a for- 



Fig. 5 Plot of log(?) versus l/T for neat resin. 
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Fig. 6 Scatter in TGA temperature data for 1% conversion. 
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changing the abscissa x and ordinate y in Fig. 5 was examined. 
Figure 7 is a plot of 1 IT versus log(^). The sample r for this 
plot also ranges from 0.91 to 0.60 for percent conversions of 
1-6%, respectively. Linear regression is again used to deter- 
mine the slopes that are inserted in the bracketed term in Eq. 
(15) to calculate E at each percent conversion. 

Table 1 lists the values of E from 1 to 6% conversion cal- 
culated with the Fortran program and both methods. Activation 
energies calculated with the standard Flynn and Wall approach 
[Eq. (14)] vary from 54 to 59 kJ/g-mol with a mean of 55 kJ / 
g-mol. These values were used to calculate the respective val- 
ues of E/RT and verify the assumptions used to derive Eq. 
(14). The observed variation in E/RT insignificantly affected 
Doyle’s constant. Using the modified Flynn and Wall method 
[Eq. (15)] the calculated activation energies vary from 66 to 
160 kJ/g-mol with a mean value of 96 kJ/g-mol . Therefore, 
using the modified method produced a much larger variation 
in E and a 75% larger mean value. A detailed comparison of 
the regression statistics revealed that both methods had iden- 
tical sample r values at the same percent conversions even 
though the actual values of E were significantly different. For 
example, at 1 % conversion the correlation coefficient for both 
methods is 0.91, but E determined with the modified method 
is 20% greater than that for the standard method. At 6% con- 
version r = 0.6 and the E calculated by the modified method 
is approximately 2.7 times larger than the value calculated with 
the standard method. However, even with a relatively low sam- 
ple correlation coefficient of 0.6 for 6% conversion, the largest 
E obtained by Eq. (14) is within 4 kJ/g-mol of the mean. Using 
the modified method the variation in E from 1 to 6% conver- 
sion is much larger relative to the mean value. However, more 
significantly, the standard errors for the modified method were 
a factor of 3—4 times smaller than the errors calculated for the 
standard method. This would indicate that the modified method 
produces a more accurate curve fit of the TGA data. The larger 
variation in E over the 1-6% conversions may also indicate 
that the modified method is more sensitive to nonlinearity of 


Table 1 Experimental activation energies 
for SC 1008 neat resin 


Percent conversion, % 

Activation energy, kJ/g-mol 

Eq. (14) 

Eq. (15) 

1 

54.0 

65.5 

2 

54.4 

67.9 

3 

53.7 

74.8 

4 

53.6 

91,4 

5 

55.6 

113.5 

6 

58.6 

159.8 

Average 1-6% 

55.0 

95.5 

Average 1, 3, and 6% 

55.5 

100.0 



Fig. 7 Plot of 1 fT versus log(^) for neat resin. 


the TGA d; ta. This nonlinearity is possibly an indication that 
the assump:ion of a first-order decomposition reaction is not 
an accurate representation for this resin. 

In practit e, the final activation energy is the average of the 
activation energies at three different percent conversions. In 
this research, 1, 3, and 6% conversions were chosen to cal- 
culate the a /erage activation energy. The final average E based 
on these va ues was 56 kJ/g-mol using the standard Flynn and 
Wall method and 100 kJ/g-mol using the modified method. 

Uncertainty in Activation Energy and Its Effect 
on Ablation Rate 

The modified method [Eq. (15)] and Eqs. (23-26) were used 
to determine the uncertainty in the calculated average activa- 
tion energy = The resultant precision Pi, bias B and total un- 
certainties Ue are, respectively, 60.0, 16.2, and 60.1 kJ/g-mol. 
The overall uncertainty of 60. 1 kJ/g-mol is 60% of the average 
E of 100 k. /g-mol. 

These results were used to evaluate the effect of uncertainty 
in activatio l energy on the ablation rate. The CMA92FLO 
code was used to determine the thermal response of a carbon - 
phenolic nozzle throat exposed to propellant gases without 
pore-pressure calculations. A baseline case was first run with 
a standard input value for E. Although the uncertainty analysis 
yielded a 60% uncertainty in E , Table 1 shows a much larger 
variation in activation energy from 1 to 6% conversion. There- 
fore, as a conservative estimate of uncertainty effects, the 
CMA code was rerun with the activation energies doubled for 
each of the three components of degrading material. The cal- 
culated ablation rates for the baseline and doubled-activation 
energy cases were compared after 30 s of a simulated solid 
rocket motor firing. Doubling the activation energy increased 
the ablation rate by 3%. Depending on several factors, a 3% 
increase in surface recession could be significant. Operating 
conditions c f the motor could cause the effect of uncertainty 
in activatio i energy to increase or decrease. In addition, 
changes in other input parameters such as enthalpy could in- 
crease or decrease the impact of activation energy uncertainty 
on ablation. 

Summary, Conclusions, and Recommendations 

Kinetic reaction parameters are used as inputs to computer 
codes that predict rocket nozzle ablation. The kinetic reaction 
parameter o activation energy is typically obtained from TGA 
data. A comprehensive literature search did not indicate that 
anyone had systematically assessed the experimental uncer- 
tainty associated with the final determination of activation en- 
ergy by one particular methodology and commercially avail- 
able instrumentation. The study described in this paper was 
conducted tc characterize the uncertainty in E and to determine 
if this uncerainty had any affect on the predicted surface re- 
cession rate of a simulated composite solid rocket nozzle. 
Carbon -phe nolic material was originally considered for this 
research; however, most data for this class of composites are 
restricted. Therefore, a neat resin, Borden SC 1008, was used. 
The uncerta nty in activation energy of a composite may be 
greater than that of a neat resin because of errors caused by 
sample noni niformity. Temperature and weight-loss data were 
gathered at leating rates of 5, 10, and 20°C/min. Activation 
energy was determined from these data. The research meth- 
odology wa; derived from established experimental uncer- 
tainty methcds, standard TGA operating procedures, experi- 
mental data, and a literature review. 

Two methods for calculating activation energy were exam- 
ined: 1) the traditional Flynn and Wall method based on the 
slope of log( q) versus 1 IT and 2) a modification of this method 
where the oidinate and abscissa are reversed. The Flynn and 
Wall approach assumes heating rate is the random variable. In 
reality the heating rate is essentially fixed and scatter is ob- 
served in the temperature data. Traditional methods of calcu- 
lating E from TGA data are actually inverse regression esti- 
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mations. Therefore, a modified method was employed in which 
activation energy is calculated with the slope determined from 
a plot of MT versus log(g). The modified method complies 
more closely with standard linear regression principles and 
produced a much larger variation in E over 1-6% conversion 
than for the standard method. Because of this variation the 
average value of activation energy using the modified method 
was 75% larger than the average value using the standard 
method. But detailed regression statistics show that the mod- 
ified method actually reduces the standard error at each percent 
conversion by a factor of 3-4 compared to the standard Flynn 
and Wall method. Therefore, the modified Flynn and Wall ap- 
proach appears to produce a more accurate curve fit of the 
TGA data and is more sensitive to data nonlinearity. 

A detailed uncertainty analysis was performed based on the 
modified method. The average activation energy for 1, 3, and 
6% conversion (£136) was calculated to be 100 kJ/g-mol with 
an uncertainty (Uf) of 60 kJ/g-mol. Based on this result, stan- 
dard activation energies for a carbon phenolic material were 
increased by a factor of 2 and used to calculate the ablation 
rate in a typical solid rocket environment. Doubling the stan- 
dard activation energy for a carbon -phenolic increased the 
predicted recession rate in a typical solid rocket environment 
by 3%. A 3% increase in surface recession could be significant 
depending on actual operating conditions of the motor and 
changes in other input parameters such as enthalpy. 

This research indicates that when linear regression is used 
to obtain information from experimental data, significantly dif- 
ferent results are obtained, depending on which variable is 
considered the dependent or random variable. This is particu- 
larly true for determination of activation energy from TGA 
data. Contrary to the approach used in the software of many 
TGA manufacturers, the activation energy should be calculated 
by a modification of the Flynn and Wall method using the 
slope determined from a plot of MT as the ordinate versus 
log(<?) as the abscissa. 
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